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The lagrangian representation of fluid motion is alternative with respect to the 
conventional eulerian description. It is focused on observing the motion of individual 
fluid particles identified by three parameters known as the lagrangian variables. The 
set of ideal-fluid lagrangian governing equations consists of four non-linear differ- 
ential equations for current positions of fluid particles as functions of the lagrangian 
variables and time [|l|]. According to Aref ([Q], p.274), ' . . . analytical tradition has 
clearly favoured the eulerian representation, and considered the lagrangian version to 
be analytically intractable'. The lagrangian analytical study of a viscous fluid is even 
more problematic because of the extreme complexity of the known expression for the 
viscous force 

This paper aims at showing that we can advance in the analytical study of the 
lagrangian fluid dynamics due to a new matrix approach proposed in [Q]. The ma- 
trix approach is based on the notion of continuous deformation of infinitesimal vec- 
tor fluid elements and treats the Jacobi matrix of derivatives of particles coordinates 
with respect to the lagrangian variables as the fundamental quantity describing fluid 
motion. The matrix ('deformation') approach is distinct from the conventional la- 
grangian ('trajectory') formulation and is, actually, similar to the tensor description 
in the elasticity theory. A closed set of governing equations is formulated in terms of 
the Jacobi matrix and studied for both ideal and viscous fluid. Recent results on the 
3-D extension of plane lagrangian solutions are reported. 

Matrix form of ideal-fluid hydrodynamic equations 

In the conventional lagrangian approach, fluid motion is described by current 
particle position X = {X, Y, Z}, considered as a function of time t and the la- 
grangian variables a = {a, 6, c}: X = X{a,t). Let us consider the Jacobi matrix 
R(a,t) = dX/da as the principal characteristic of fluid motion. The evolution of 
an infinitesimal material vector AX corresponding to an increment of the lagrangian 
variables da is then given by dX = Rda. For an ideal fluid, the Jacobi matrix is 
shown to satisfy the following set of equations [Q] : 

det R = det Rq, R("R — R'^R* = S, dRnm/dak = dRnk/dam, (1) 

where Rq = R(a,0), R^ is the transposed Jacobi matrix, subscript 't' denoting 
the time derivative, S = [cijkS^) is a time-independent antisymmetric matrix of the 
Cauchy invariants, {a,b,c} = {01,02,03} for the indexed notation. The Cauchy 
invariants are the integrals of the conventional 3-D lagrangian equations related 
to the initial vorticity at t = by 5^(a) = {f2o, Vcfc) det Rq. Set ([I]) comprises 
the matrix analogues of the lagrangian continuity and vorticity equations along with 



the consistency condition. This condition allows us to reconstruct the particle tra- 
jectories X{a,t) by a solution of ([l]) for R uniquely through integration over the 
lagrangian variables: 

X{a,t)=f^^R{q,t)dq. (2) 

The specific features of set (|l|) consist in the following. The continuity equation 
takes the form of an algebraic constraint on the Jacobi matrix. The equation of mo- 
tion is formulated as a nonlinear matrix differential equation in time only, where the 
derivatives with respect to the lagrangian variables do not appear. The accompanying 
consistency condition is linear and does not include the time derivatives. Altogether, 
set (|I|) admits a flexible formulation of the problem of interest and allows one to take 
advantage of the powerful machinery of matrix calculus. 

The expression for vorticity takes the form 

/2 = (det Ro ) " ^ 5* dX/dai (3) 

and represents decomposition of vorticity into components along the reference vec- 
tors dX /dtti of the frozen-in 'liquid' coordinate frame of the lagrangian variables. 
Accordingly, the meaning of the invariants is that they are proportional to the 
contravariant components of vorticity in the 'lagrangian' reference frame. As seen 
from the formula for its Cartesian components i7„ = detRg S\ the evolution 
of vorticity is described by the Jacobi matrix directly. 

Thus, by virtue of the facts that set (|I|) is closed and one can reconstruct any 
required characteristic of the flow by its solution for R it follows that the matrix 
approach based on set (|l]) provides a complete description of motion of an ideal fluid. 

A method for constructing solutions to the set of matrix equations which corre- 
spond to 3-D rotational flows with precessing vorticity is developed in [Q]. It allowed 
us to find and investigate a new class of 3-D exact matrix solutions that would be ex- 
tremely difficult to derive from conventional formulations. The solutions investigated 
in depend not only on a set of parameters but on several arbitrary functions of the 
Lagrangian variables and incorporate the known 2-D solutions such as the Ptolemaic 
vortices [^], the Gerstner waves and Kirchhoff's vortex as particular cases. Below 
we report results of an alternative technique to construct 3-D solutions of the matrix 
equations (|l|). 

3-D shearing and stretching of plane motion 

The matrix formulation provides a compact 'block' representation of 3-D lagran- 
gian equations and allows us to study new 3-D extensions of known 2-D rotational 
motions. Suppose that X{a,b,t),Y{a,b,t) represent a solution of 2-D lagrangian 
hydrodynamic equations. The corresponding Jacobi matrix takes the form 
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1 



(4) 



It can be proved that along with matrix set (|I|) is satisfied by the transformed 
matrix 
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(5) 



which includes an arbitrary function of the lagrangian variables H{a, b) and an arbi- 
trary function of time f{t) ^ 0. The particle trajectories obtained from @ using @ 
take the form 

X {a, b, J fdt) Y {a, b, J fdt) ,w,-2H,^l ,M 

^ = fTj2 ' ^ = JTj2 ' ^ = /[^(«'^)J/ dt + cj (6) 

and represent a 3-D motion with the Z-component of velocity 

V, = [H{a,b) + fZ]/f. (7) 
The Cauchy invariants of the 3-D solution are given by 

S'^ = Hh, S"^ = —Ha, = XtaXh — XfhXa + YfaYf, — Yt^Ya, (8) 

where 5''^ remains the same as for the 2-D solution (Q), but both and are now 
non-zero. The vorticity corresponding to @ is obtained from (||) using (|8]), it 
becomes 3-D and time-dependent: 

n = {f-'dH/dY, -f-'dH/dX, S'f/iXaY, - X^Ya)] . (9) 

The transformation according to (||), allows one to obtain non-trivial exten- 
sions of any given 2-D solution. As can be seen from (0), (^), the contributions of 
arbitrary functions H{a,b) and f{t) to the resulting 3-D motion are a shear flow 
along the Z-direction and stretching of the Z-component of vorticity, respectively. 
We shall illustrate the effect of these two factors separately using the 2-D Ptolemaic 
vortices [Q, ^ as the initial solution. The lagrangian expression for the flow inside 
the 2-D Ptolemaic vortex is given in complex variables [Q] : 

= 2-1/2 (X + iy) =eexp(-2iwt) + F(0 f^j. |^| < (iq) 

where W is the complex coordinate, ^ = 2^^^^ {a + ib) is the complex lagrangian 
variable, ^ denotes the complex conjugate variable, function F{^) is analytical within 
1^1 < 1 and satisfies there the condition < 1. 

3-D Ptolemaic vortices with swirl 

In order to focus on the effect of the Z-shear, we apply transformation (Q) to ( [T0| ) 
and eliminate stretching by setting f{t) = 1. The resulting flow is most expediently 
described in complex form similar to the initial solution ([T0|): 

H^ = eexp(-2ia;t)+F(e), Z = H{^,Ot + c. (11) 



Here, H{^, ^) remains real-valued. The particle trajectories are regular circular he- 
lices. Although the lagrangian expression for the Z-component of the flow is ele- 
mentary, the eulerian velocity for cannot be written explicitly and is evidently 
non-trivial. The Cartesian components of vorticity are derived from @, (|ll|): 
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(12) 



To have a geometrical idea of the structure of vortex lines, consider their equation in 
the space of the lagrangian variables da/S^ = db/S'^ = dc/S^. The substitution 
of the Cauchy invariants from (^ yields da/Hh = —db/Ha = — dc/(4a;), then 
H = const along the vortex lines. If the contours of H{a, b) = const are closed 
inside the vortex (for \a + ib\ < 1), the vortex lines are windings on the cylindrical 
surfaces in the space of a, b, c, for which these contours serve as directing curves. 
The shapes of vortex lines in Cartesian space are topologically similar, although they 



experience periodic distortions, since the smooth mapping defined by (11) is time- 
dependent. Figure |I] shows the material (moving with fluid particles) vortex lines of 
flow (Q) for I I = 0.9, and the surface of the vortex tube for |^| = 1: (a) at t = 0, (b) 
at t = 2tt- here F = a{l + /?)^/^ a = 0.11, /? = 1.2, H 

Such an important solution of 2-D lagrangian equa- 
tions as the Gerstner waves |[T], ^, can be extended by 
analogy with (11): 



W = C + Aexp[i{kC + crt)], Z 



(13) 



Noteworthy is the fact that (pj|) satisfies not only the 
governing equations but also the boundary conditions 
at the free surface Im ^ = 0. The resulting 3-D flows 
have helical particle trajectories and generally curvi- 
linear vortex lines, their shape depending on the real- 
valued function H{S,, ^). 

Stretched Ptolemaic vortices 
Assume that the Z-shear in (|5|)-(^ vanishes with 
H = 0, and consider the effect of the arbitrary function 




Figure 1 : Vortex lines (solid) 
and vortex tubes (dotted) for 
dlTI ) (see the above text for 
detail). 



f{t) proceeding from 2-D Ptolemaic vortices (10). We 
shall demonstrate that the extension according to (^) 
leads to stretching or contraction of the vortex by a non- 
stationary axisymmetric strain flow. The particle trajectories obtained from (|l^) for 
H = 0, f{t) / take the form 



W = [e exp {-2iu: ffdt)+ F{0] , Z = f{t)c. 



(14) 



For monotonic f{t), all the particles trace convergent / divergent spiral paths. Solu- 
tion (14) has one non-zero Cauchy invariant 5^ = —Auj and the only component 



of vorticity Qz = — (l — |-F'P) which depends on the lagrangian variables 
and time. It is the time variation of vorticity in proportion to f{t) that reflects the 
effect of concentration and intensification of vorticity (for d|/|/di > 0) due to the 
stretching of vortex tubes along the Z-direction. 

The manifestation of the strain flow becomes evident for the isolated vortex (jlj) 
streamlined by an irrotational fluid. We now assume that the 3-D rotational motion 
( [T4| ) is concentrated within a vortex core, e.g. within the vortex tube |^| ^ 1. It proves 
possible due to a specially developed technique of potential continuation to construct 
an exterior 3-D potential flow which agrees with the rotational motion across the 
boundary of the vortex. We formulate a system of parametric matrix equations for 
the exterior potential velocity, provided that parameterization of the coordinates in 
the outer region is given. The set of parameters for the outer potential region is most 
expediently chosen in accordance with the parameterization of the interior by the 
lagrangian variables ([l^. Using the matrix equations of potentiality we obtain the 
potential velocity paiametrically at the point of the outer region 

W = r^'^ [i exp {-2iu f fdt) + F{l/0] , Z = f{t)c for |C| > 1 (15) 

as a function of parameters ^, c : 

V = 2-^l\V^ + iVy) = -2i../V2e-2-//d7^_ 1 VF///, Vz = fc, (16) 



where / stands for the time derivative of /. All the components of ( |16| ) agree with 
the velocity obtained from (jlj) across the vortex boundary |^| = 1. Unlike the 
vortex core, for the outer region ^, c are not the lagrangian variables of the poten- 
tial flow and play the role of formal parameters in ([Tsl), ([l6| ). As the distance from 
the vortex increases (|^| oo), the potential velocity (|16D tends asymptotically to 
the non-stationary axisymmetric strain flow, the intensity of which is proportional 
to /: Vr >i — r//(2/) + 0(l/r), = Zf/f, where r is the cylindrical radius. 
It becomes clear that the stretching or contraction of the vortex is influenced by 
the external strain and is essentially kinematical. Figu re ^ shows the boundary of 
the material section |^| = 1, |c| ^ 4 of the vortex (|T^([I^ for f{t) = t and 
Fil) = Ci)"^ + - C2y^], a = 0.83, Ci = 2.2, Ci = 2.09 + 0.68i. 

We have demonstrated the effects of 3-D shear and axisymmetric strain on 2-D 
lagrangian solutions in succession. As seen from general expressions (^, (^, the 
extension of 2-D solutions allowing for their concurrent contribution is also simple 
and results in a variety of interesting 3-D flows. 

Lagrangian equations for a viscous fluid 

The notion of the frozen-in lagrangian coordinates and the matrix formulation 
are useful not only for an ideal fluid but also for a viscous fluid. We have seen that 
the constancy of the Cauchy invariants plays an important role for an ideal fluid. 
It is natural to question how they change due to viscous force. We shall derive an 
equation for which serves as the lagrangian analogue of the Helmholtz equation 
for a viscous fluid. 



The lagrangian equations deduced by a direct 
change to the lagrangian variables in the Navier- 
Stokes equation are known from the literature |j3|]. 
However, their structure is so complicated that no 
progress in their analytical study or applications 
has been achieved until now. 

For an alternative approach, we consider the 
column vector S = {S^,S'^,S^} composed of 
the Cauchy invariants. It is related to the velocity 
vector by 5 = vota{R^Xt), where (rotaA)j = 
CijkdAj /dak- The change of S occurs due to the 
viscous force exclusively and is governed by the 
diffusion-type equation 

dS/dt = -7] rota [g rota (g S/D) /D] , (17) 
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Figure 2: Stretching of a material 
section of the vortex (|l3)-(|l6l) 
in a strain flow. 



where t] is the kinematic viscosity, g = {go.) = R^R is the metric tensor reflect- 
ing the distortion of the frozen-in 'liquid' coordinate system, D = (detg)^/^. This 
equation complements set (|I]) of the ideal-fluid matrix equations. Together, they con- 
stitute a closed system of equations that fully describes the motion of a viscous fluid, 
provided that 5* in (|T]) are considered as functions of time. 
For a 2-D motion = S'^ = and (|l^) reduces to 
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(18) 



where 5 = 5''^. Equation (18) can be simplified asymptotically for low or high 
viscosity. For the case of low viscosity (Re <C 1), the tensor components gik can 
be approximately replaced by their time average obtained from the inviscid problem. 
For high viscosity (Re S> 1), the term dS/dt on the left-hand side o f (|I^ ) can asymp- 
totically be neglected, so that time plays the role of a parameter in (|18[) (note that we 
should nevertheless keep the time derivatives in ([T])). More profound study of (|T7|), 



(18) and their applications are beyond the scope of this report. 
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